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We study chemical equilibration of quarks and gluons in central nuclear collisions at 
RHIC and LHC energies. The initial quark and gluon densities are taken from earlier 
studies as well as from recent perturbative QCD estimates and are then evolved via 
rate equations coupled to longitudinally boost-invariant fluid dynamics. We find that, 
for RHIC initial conditions, the lifetime of quark- gluon matter is too short in order for 
the quark and gluon number densities to chemically equilibrate prior to hadronization. 
In contrast, at LHC energies chemical equilibration is complete before the system 
hadronizes. Entropy production due to chemical equilibration can be as large as 30%. 

x. 

1 Introduction 

Nucleus-nucleus collisions at ultrarelativistic energies probe the properties of nuclear matter 
under extreme conditions Q. Of particular interest is the question whether ordinary nu- 
clear matter undergoes a (phase) transition to quark-gluon matter, as predicted by lattice 
calculations of quantum chromodynamics (QCD) ||. 

A nuclear collision can be viewed as a sequence of nucleon-nucleon collisions. At suffi- 
ciently high energies, multi-particle production leads to the formation of a region of high 
energy and particle number density. With increasing beam energy, multi-particle production 
through processes where the partons inside the nucleons directly interact with each other, 
becomes more and more important. The dominant partonic particle-production mechanism 
is so-called mini-jet production 0. 

If the momentum transfer in these interactions is sufficiently large, mini-jet production is 
reliably computable within perturbative QCD (pQCD). It was estimated [[| that about half 
of the total transverse energy created in Au + Au collisions at the Relativistic Heavy-Ion 
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Collider (RHIC) at BNL resides in mini-jets, while at CERN's Large Hadron Collider (LHC) 
the transverse energy created is almost exclusively due to mini-jets. 

The transverse energy can be used to estimate the energy densities created in the initial 
stage of the collision. If matter is in thermodynamical equilibrium, both for RHIC and LHC 
energies the corresponding energy densities are found to be large enough for nuclear matter 
to be in the quark-gluon phase. The remaining question is whether high-energy density 
matter formed in ultrarelativistic nuclear collisions lives sufficiently long enough to actually 
reach thermodynamical equilibrium, so that it can be identified with the quark-gluon phase 
seen in lattice QCD calculations. 

Thermodynamical equilibrium means that matter is in thermal, mechanical, and chemical 
equilibrium. In general, the question of thermodynamical equilibration can only be decided 
with microscopic transport models ||. In this paper, we study a simpler problem which 
allows us to use a macroscopic transport model. We assume that matter reaches thermal and 
mechanical equilibrium after a proper time tq ~ 0.1 — 0.25 fm/c. This is justified given the 
fact that the rate for elastic collisions between quarks and gluons, which establish thermal and 
mechanical equilibrium, is much larger than for inelastic collisions which establish chemical 
equilibrium. Thermal and mechanical equilibrium will be referred to as kinetic equilibrium 
in the following. 

Under these assumptions, and given initial values for the energy density and the quark, 
antiquark and gluon number densities, we can then employ ideal fluid dynamics to study 
the subsequent evolution of the kinetically equilibrated quark-gluon phase, coupled to rate 
equations which determine the chemical composition of the system away from chemical 
equilibrium. This problem has been previously studied by Biro et al. 0, and by Srivastava 
et al. 0, |S|. As in these previous studies, we assume that chemical equilibration is driven 
mainly by the two-body reaction gg <-> qq, and by gluon multiplication as well as fusion, 
99 999- We also terminate the evolution when the plasma reaches the hadronization 
energy density, = 1.45 GeV/fm 3 . For a completely equilibrated plasma, this corresponds 
to a hadronization temperature of T h ~ 0.17 GeV. We do not study the evolution in the 
mixed and purely hadronic phases. 

Our motivation to reinvestigate this subject is the following. First of all, the question 
whether the quark-gluon phase created in ultrarelativistic heavy-ion collisions is chemically 
equilibrated is highly relevant for the experiments commencing at the RHIC collider in the 
fall of this year. Since |5], |7j, were published, newer results on mini-jet production in 
nuclear collisions at RHIC and LHC energies have been obtained [|J. We therefore decided 
to compare the evolution computed with initial conditions obtained from the so-called "self- 
screened parton cascade" (SSPC) model (used in [0, §|) with that employing the more recent 
pQCD estimates of Q. 

Second, the authors of || [7], |[ simplified the rate equations using an approximate, so- 
called "factorized" phase-space distribution function for the quarks and gluons (see below). 
In our treatment we use the full distribution function in the rate equations, and so are in a 
position to assess the validity of the factorization assumption. 

Third, baryon stopping is non- negligible in nuclear collisions at SPS energies ||, and there 
is mounting evidence from theoretical studies that, even at RHIC energies, the midrapidity 
region is not completely net-baryon free ||10|| . Our analysis therefore also accounts for non- 
zero net-baryon number. 
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Finally, the numerical algorithms used here are different from the ones employed before 
]7]. § . Thus, they constitute an independent check on the validity of the conclusions reached 
previously. 

This paper is organized as follows. In section 2 we discuss the macroscopic transport 
equations. Starting from a single-particle phase-space distribution in kinetic equilibrium, 
we show that these equations are given by ideal fluid dynamics, coupled to rate equations 
which determine the densities of the individual particle species. We also prove that the 
entropy never decreases during the time evolution of the system. In section 3 we discuss the 
initial conditions for the transport equations. Section 4 is devoted to purely longitudinal 
boost-invariant expansion. We show that the difference in the equilibration process is small 
when using the full phase-space distribution function instead of the factorized distribution 
function of [|, |7|, ||. We compare results obtained with the initial conditions of the SSPC 
model J7], |5J with those when using pQCD estimates M. We also study the sensitivity of the 
equilibration process on the initial time, tq, and the strong coupling constant, a s . In section 
5 we include (cylindrically symmetric) transverse expansion as well, and present results for 
the SSPC model. We conclude in section 6 with a summary of our results and an outlook. 
Our units are Ti = c = k B = 1, and the metric tensor is g^ v = diag (+, — , — , — ). 

2 Macroscopic transport equations 

In this section, we discuss the macroscopic transport equations on which our results are based. 
We assume that elastic collisions between quarks (antiquarks) and gluons are sufficiently 
frequent to establish kinetic equilibrium. Ideal fluid dynamics can then be invoked to follow 
the evolution of the energy and momentum densities in the system, while rate equations 
for the inelastic reactions gg <-> ggg and gg <-> qq determine the number density of quarks, 
antiquarks, and gluons. 

2.1 Phase-space distribution 

Once kinetic equilibrium is achieved, all particles in an (infinitesimal) volume element at 
space-time point x = x^ = (t, x) have the same temperature, T(x), and move with a common 
average 4- velocity, u{x) = u^(x). In that case, the single-particle phase-space distribution 
for particles of species i assumes the same form as in (local) thermodynamical equilibrium, 
except that the fugacity (or the chemical potential) is not equal to its equilibrium value: 



where k = k^ = (Ei, k) is the 4- momentum; Ei = y k 2 + mj is the on-shell energy of particles 
of species i with 3-momentum k and rest mass m.j. The fugacity Xi(x) = exp[p,i(x)/T(x)] 
controls the number density of particle species i, //j is the chemical potential. 9i = ±1 for 
fermions or bosons, respectively, dj denotes the number of internal degrees of freedom for 
particles of species i (like spin, isospin, color, etc.). 



fi(x, k-u) = Xi(x) 



di 1 




(2tt) 3 exp [k ■ u(x)/T(x)} + X^x) 6, t ' 
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In [||, [j], §], instead of (|T|) a so-called factorized distribution function was used, in which 
the fugacity Aj in the denominator is approximated by 1, 

ft c (x, k-u) = \(x) j^- , Lv^a • ( 2 ) 

(27r) ,:i exp [k ■ u{x)f 1 (x)\ + Ui 

We will comment on the validity of this approximation below. 
2.2 Energy- momentum conservation 

With a distribution function of the type ([I]), the energy-momentum tensor assumes the 
so-called ideal-fluid form 

T^(x) = V / — k" k u fi(x, k-u) = [e(x) + p(x)\ u^x) u u (x) - p(x) , (3) 
i J Ei 

where 
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e(a;) = W tfkEififaEi) , p(a?) = £ / d 3 k^/,(x,^) (4) 

i i ' 

are the energy density and pressure in the local rest frame of a fluid element moving with 
4- velocity u M . 

While it is clear that thermal equilibrium requires all particle species to have the same 
temperature, it is less obvious why the phase-space distribution ([I]) ensures mechanical 
equilibrium as well. To see this, consider the tensor decomposition T^ v as given by (|3|). In 
the rest frame of the fluid element, the pressure is completely isotropic, p = T u /3 = T xx = 
T yy = T zz , which is synonymous to mechanical equilibrium. If ([J) depends on more than 
one 4- vector u, additional tensors would appear on the right-hand side of (^), such as u^u\, 
and the diagonal components of the stress tensor T 4 - 7 would no longer be identical. Hence, 
the common 4-velocity u in ([!]) ensures mechanical equilibrium. 

In our case, matter consists of massless gluons, quarks, and antiquarks, i = g, q, q, with 
energy density 

e = e g + e q + eg , (5) 

A poo ^3 

rj POO */3 

^ = KT'^-J dz^—, (7) 
2n z Jo e z + X q 

e q = XgT'A Hdz ' 



2tt 2 Jo e z + X q 



where d g = 2(N% — 1) is the number of internal degrees of freedom for gluons, N c = 3 is the 
number of colors, and d q = 2N c Nf, is the number of internal degrees of freedom for quarks 
and antiquarks, with Nf being the number of massless flavors. Throughout our analysis, 
we use Nf = 2.5, mimicking the effect of the nonzero mass of the strange quark by taking 
d s = 0.5 d q . 
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This approximation interpolates between the region of high temperature, T 3> m s ~ 150 
MeV, where Nf ~ 3, and the region of low temperature, m u ^ < T < m s , where Nf ~ 2. 
Since we follow the time evolution of the temperature from very high T down to T ~ T c ~ m s , 
taking Nf = 2.5 is not a particularly good approximation. We nevertheless make this choice 
to be able to compare our results to earlier work |], [7], ||. Also, the chemical reaction rates, 
cf. eqs. fl22T) below, which appear on the right-hand side of the rate equations fll9l) - fl2l|) , are 



more complicated in the case of massive particles ||11|| . The extension of our present study 
to nonzero strange quark mass is, however, important and will be pursued in a subsequent 
publication JTJ| . 



For massless bosons, the equilibrium value of the fugacity, Af, is always equal to 1. The 
reason is the following. For bosons, the equilibrium value of the chemical potential, /4 q , 
has to be smaller than the rest mass, /if 1 < Wj. On the other hand, bosons which carry a 
conserved charge, /if ^ 0, always come in pairs with their own antiparticles, e.g., tt + and 
7r~, K and K, etc.. Without loss of generality, we can therefore choose /if > (with the 
chemical potential of the associated antiparticle /if = —/if < 0). Then, for massless bosons, 
/if = 0, such that Af = 1. 

For massless fermions, there is no such constraint, although Af = 1/Af in all cases and 
Af = Af = 1 for vanishing net-baryon density. 

The pressure is related to the energy density through 

V=\, (9) 

the well-known equation of state for an ultrarelativistic ideal gas. Note that, in order to 
derive this relation, we only required the system to be in kinetic equilibrium, i.e., eq. (0) is 
valid even when the system is not in chemical equilibrium! 

With the factorized distribution function (Q), the integrals in (^[) - ([8]) can be performed 
analytically, as they no longer depend on the fugacities. The energy densities simplify to 

4 ac = a 2 A 9 T 4 , el ac = b 2 X g T 4 , e<f c = b 2 A 9 - T 4 , (10) 

where a 2 = 8n 2 /15 and b 2 = 7n 2 Nf/4:0. The equation of state (|9]) remains valid. 

In Fig. |l] we show the ratio e^ ac /ej for quarks and gluons as a function of the fugacity 
Aj. This ratio does not depend on the temperature. For bosons, the factorized expression 
tends to overestimate the correct result, for fermions, it underestimates it. In the interval 
< Aj < 1, the error is maximized at Aj = 0, and decreases monotonically as Aj — > 1. At 
Aj = 0, it is 8% for bosons, and 5% for fermions. 

Energy and momentum is locally conserved, 

d^(x) = . (11) 

With (0), these are the equations of ideal fluid dynamics, which can be used to compute the 
evolution of the energy density and the fluid 4-velocity of the system. The fluid-dynamical 



equations (JTT|) are closed by specifying the equation of state of matter under consideration, 
i.e., the pressure as a function of energy density, p = p(e). In our case, this equation of state 
is simple, see eq. @. 
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Figure 1: The ratio ef^/ej for gluons (solid line) and quarks (dashed line) as a function of 
the fugacity Aj. 

2.3 Rate equations 

The 4-current of the number of particles of species % is 



is the number density of particle species i in the rest frame of a fluid element. This density 
is controlled by the value of the fugacity, Aj. In general, the fugacities for the different 
particle species are independent thermodynamic variables, and one has to specify additional 
equations of motion which determine the number densities. These are the rate equations. To 
close the coupled set of fluid-dynamical equations and rate equations, one has to be specify 
an equation of state which, in general, also depends on the number densities of the various 
particle species, p = p(e,n\,n2, • ••)■ 

The reason why there is actually no such dependence for the system considered here, eq. 
is that the quarks and gluons are considered to be massless. This has the important 
consequence that, in our case, chemical non-equilibrium does not affect the evolution of 
the energy and momentum densities of the fluid! This certainly changes in the case of a 
nonzero strange quark mass, although, for large temperatures T 3> m s , the dependence of 
the pressure p on n s is relatively weak. 

The number densities of gluons, quarks, and antiquarks in the rest frame of a fluid element 

are 




(12) 



where 




(13) 




(14) 



(15) 
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Figure 2: The ratio nf^/rii for gluons (solid line) and quarks (dashed line) as a function of 
the fugacity Aj. 



A„-T 3 dq 



2tt 2 



dz 



e z + Xn 



With the factorized distribution function (Q) these expressions simplify to 



rig" = a 1 X g T 



n 



fac 



n 



(16) 



(17) 



where a x = 16C(3)/vr 2 and b x = 9 C(3) N f /(2n 2 ). 

In Fig. ^, we show the ratio n l ^ c /rii as a function of fugacity. As for the ratio of energy 
densities, Fig. [I], this ratio does not depend on the temperature, and the factorized expression 
overestimates the correct result for bosons, while it underestimates it for fermions. Again, 
the error decreases as a function of Aj, but it is about twice as large as for the energy 
densities. The maximum error (at \ = 0) is about 20% for gluons and 10% for quarks. 

In the absence of chemical equilibrium, the number of particles of species i is determined 
from the rate equation 

d^N?(x) = Ri(x) . (18) 

In chemical equilibrium, Ri vanishes. In our case of a quark-gluon system, we assume that 
the Ri are determined by the reactions gg <-> ggg and gg <-> qq. The rate equations (|T8|) can 
then be written in the form 101 



d,{n q u») 
d^nqu 11 ) 



R 3 n s 



R 2 n q 



R 2 n Q 



2R 2 n„ 



n q rig n 2 g * 
n q n q n 2 gi 
n q n q n 2 g ' 
h q h q n 2 g/ 



n q n q n 2 g 
n q n q n 2 g 



(19) 
(20) 
(21) 
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where the hi are the number densities computed at the same temperature T, but for Aj = 1. 
The terms R2 and R3 are given by M 



R 2 ~0.2AN f a 2 s X g T\n(1.65/a s X g ) , R 3 ~ 2.1 a 2 s T^2X g - X 2 g , (22) 

with a s being the strong coupling constant. 

We take a s = 0.3 throughout this analysis, unless otherwise stated. Note that, while 



the densities on the left-hand side of eqs. (flTf ) - (|21|) in general contain the full phase- 
space distribution function ([!]), the right-hand sides of these equations and the terms (|22| ) 
are derived assuming the f actor ized distribution (H) ||. This is in principle inconsistent. 
At least close to equilibrium, however, the right-hand sides of the rate equations (|T9| ) - 
([HP are small, and we are allowed to neglect this inconsistency. On the other hand, far 
away from equilibrium we expect other contributions to the chemical reaction rates to be 
more important than the difference between factorized and full distribution functions. Such 
contributions are, for instance, multi-gluon processes gg — > ng [Tj|. In this work, we only 
include the lowest-order rates (|2T 



In chemical equilibrium, the right-hand sides of the rate equations ( |l9|) - (21) vanish by 
definition, and we obtain 

W/) = 0, i = 9,q,q- (23) 

We can now derive a relation for the time evolution of the ratio of the parton number density 
rii to its corresponding equilibrium value nf. By writing rii = nfni/nf on the left-hand 
sides of the rate equations (0) - (|2lD , we derive with eq. (|23"D the condition 




= 4 • (24) 

This means that, for R^jnf > 0, the comoving time derivative of the ratio of density 
to equilibrium density is positive, and thus this ratio is bound to grow with increasing 
proper time in each fluid element. Vice versa, this ratio decreases with time if Ri/nf < 0. 
In other words, in the rest frame of a fluid element, ni/nf 1 cannot decrease (increase) if 

iy< q >o(<o). 

2.4 Entropy 

Local thermodynamical equilibrium implies the conservation of entropy. If the system is not 
chemically equilibrated, the entropy is bound to grow, until chemical equilibrium is reached 
|7], |H| . To see this, we make the following observation: the form of the distribution function 
(H|) or (|2|) is exactly that of a distribution function in local thermodynamical equilibrium, 
except that the fugacities, Aj, assume values which are not the same as in chemical equilib- 
rium, Af 1 . Note that in this case there are no further dissipative terms |15|] in the tensor 



decomposition of (§) and (12) [|16|]. Nevertheless, the system is not in thermodynamical 



equilibrium, and the actual values of e, p, and rij will differ from the equilibrium values, 
e eq , p eq , and nf 1 . 

Since the distribution function has the same form as in thermodynamical equilibrium, we 
can imagine each fluid element to actually be in chemical equilibrium with a (local) particle 



bath which determines the value of the fugacities instead of the rate equations (|Tq) . In 
other words, we can imagine local changes of particle number to be induced by a change of 
the parameters of the particle bath instead by chemical reactions within the fluid element 
itself. As long as these changes happen faster than those of the macroscopic fluid variables, 
one can view the fluid element to remain in (local) thermodynamic equilibrium, and apply 
thermodynamical relationships. 

The assumption that microscopic reaction rates are much larger than the rate of change 
of macroscopic fluid variables is certainly an over-idealization in view of the characteristic 
time scales in relativistic heavy- ion collisions. So is, in fact, the ideal fluid limit which 
assumes complete (local) thermodynamical equilibrium. Nevertheless, our point of view is 
that ideal fluid dynamics still offers valuable information about the collective behavior of the 
system. Our approach is the simplest possible extension of ideal fluid dynamics. A better 
approximation would be to solve the equations for dissipative relativistic fluid dynamics in 



the presence of chemical reactions [ lTfl . In general, chemical non-equilibrium gives rise to 
additional diffusion terms which do not appear in our treatment. So far, the complexity 
of these equations has discouraged attempts to apply them to the description of nuclear 
collision dynamics. 

Accepting these caveats, we now proceed by contracting (|TT|) with u v , 

= u u d^ v = u-de + (t + p)d-u , (25) 

and use the first law of thermodynamics, 

de = T ds + ^ drii , (26) 

i 

(s is the entropy density), as well as the fundamental relation of thermodynamics, 

e + p = T s + ^2fiirii , (27) 

i 

to derive 

= T (u ■ ds + s d ■ u) + J2 ■ drii + riid ■ u) . (28) 



Equations ( |18"D and (|12"D, together with the definition of Aj, can be used to conclude that 

<V = "E lnA ^ , (29) 

i 

where = s/. In chemical equilibrium, Ri vanishes, and entropy is conserved. 

Consider now the case of gluons, or quarks when the net-baryon number is zero. In this 
case, the equilibrium value of the fugacity Af 1 = 1, i.e., /zf 1 = 0. Now assume that the 
actual particle number density is smaller than the equilibrium value, corresponding to 
Aj < A| q = 1. The rate equations drive Hi towards equilibrium, i.e., the right-hand side has 
to be positive, Ri > 0, in order to produce particles of species i. Since In Aj < 0, the change 
in entropy is positive. On the other hand, if is larger than in equilibrium, Aj > Af 1 = 1, 
the rate equations reduce the number of particles of species i, i.e., the right-hand side is 
negative, Ri < 0. Again, entropy increases, since lnA« > 0. 
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Now suppose that Af 1 > 1, for instance for quarks when the net-baryon number is 
positive. Naively, one would think that the previous argument fails in this case, and one 
might worry whether this could lead to a situation where entropy actually decreases towards 
its equilibrium value, which contradicts the second law of thermodyamics. This is, however, 
not the case. In a relativistic system, one can create quarks only in pairs with antiquarks. 
In order to conserve the net-baryon number, d^N^ — Nq) = 0, the right-hand side of the 
rate equation for the quark number density, R q , must be equal to the right-hand side of the 
rate equation for the antiquark number density, Rg, see eqs. fl33]) - Cj34|) below. 

In eq. (^), these two terms in the sum over i can then be combined to — (fi q + fiq)R q /T. 
In equilibrium, /i® q = — /i| q , such that fig + fig = 5fig + 5fig, where = /Xj — /i? q . If the quark 
(and consequently, the antiquark) number density is smaller than its equilibrium value, 
5fi q + 5 fig < 0, the right-hand side of the rate equation is positive, R q > 0. Consequently, 
the sum of the quark and antiquark term on the right-hand side of (|2D| ) is positive, leading 
again to an increase in entropy. The argument is similar in the case when the quark number 
density is larger than in equilibrium, or when Af 3 < 1. In all cases, the entropy increases. 



3 Initial conditions 

In order to solve the macroscopic transport equations discussed in the last section, one has 
to specify the initial conditions. In this paper, we use results from two different approaches, 
both for RHIC and LHC energies. 

The first approach is the SSPC model employed in 0, |§. We decided to use these 
initial conditions for two reasons. First, the values for initial energy and parton densities 
obtained in this approach constitute an upper bound of what is expected to be created in 
ultrarelativistic nuclear collisions at RHIC energies. Second, it allows us to directly compare 
the results of |7|, |[ with ours and point out possible discrepancies. 

The second approach is that of ||. This approach is purely based on pQCD, and thus 
does not contain the contribution from the soft background, which constitutes about half of 
the total transverse energy at RHIC energies. In this sense, at RHIC this second approach 
constitutes a lower bound for the initial energy and parton number densities. 

Table || contains the initial values for energy and parton number density for each approach 
and collision energy. The SSPC values are taken from pj, while the values shown for the 
pQCD approach |4| are computed from the dN/dr] and dE T /di] values near midrapidity 
(—0.5 < 7] < 0.5) using a volume element AV^ = 7rR 2 T Ari ||18|| , such that 

d^T„ dV dN,_ Q dV . _ 

—— = e — = e T irR , -— = Ui — = n^nR , i = g,q,q, (30) 
dr] dr/ dr/ dr] 

where R = 1.12 A 1 / 3 fm. From the four values eo, n°, and n^, using eqs. (|5|) and ( |T^ ) 
(P~6|) one can unambiguously extract T , Ajj, A°, and A^, and, from the latter, the values 



n° / n° given in the last three columns of Table [l]. 

The initial time r in the SSPC and pQCD approaches is chosen as r = 1/po, where p is 
the infrared momentum cutoff required to regularize the mini-jet production cross sections. 
The values for tq given in Table [I] arise from the cutoffs po ~ 0.8 GeV for the SSPC approach 
and p ~ 2 GeV for the pQCD approach. The transverse momentum p T > p of a produced 
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Table 1: Initial conditions for the fluid-dynamical expansion phase in central Au + Au 
collisions at BNL RHIC and Pb + Pb collisions at CERN LHC energies from the SSPC 
model and pQCD. 



Approach 


Energy 


To 




T 


n° g /h° g 


n° q /n° q 


n°g/n°g 






(Im/c) 


[\je V/tm ) 


(be V ) 




SSPC 


RHIC 


0.25 


61.4 


0.66 


0.34 


0.064 


0.064 


SSPC 


LHC 


0.25 


425 


1.01 


0.43 


0.082 


0.082 


pQCD 


RHIC 


0.10 


23.9 


0.889 


0.042 


0.0077 


0.0048 


pQCD 


LHC 


0.10 


1057 


1.09 


0.99 


0.064 


0.061 



mini-jet determines the time scale for the respective parton to come on mass-shell. Only 
after the partons are on their respective mass-shell, elastic scattering processes between 
them will drive the system towards kinetic equilibrium. How the system approaches kinetic 
equilibrium has to be studied in the framework of kinetic theory; our approach does not make 
any statements about this pre-equilibrium stage of a heavy-ion collision. In the following, we 
take To = 1/po, which is certainly the earliest possible time for the system to reach kinetic 
equilibrium. In the next section, we shall also use larger (and thus more conservative) values 
for ro, in order to test the sensitivity of our results on this important parameter. 

For RHIC initial conditions, the energy density in pQCD is much smaller than for the 
SSPC model, despite the smaller initial kinetic equilibrium time r . The reason is the absence 
of the soft background contribution in the former approach. For LHC, pQCD gives larger 
values for the initial energy density than the SSPC model. If scaled to the same initial 
proper time, however, the initial energy densities are approximately equal. This reflects the 
fact that the soft background contribution is rather small at LHC energies. Nevertheless, 
despite the similarity of the initial energy densities, in pQCD the gluons are close to chemical 
equilibrium at LHC energies, while this is not the case for the SSPC approach. Note that 
there is a small initial net-baryon number density (n° — n^)/3 in the pQCD approach, as 
opposed to the SSPC model where vanishing net-baryon number is assumed. 



4 Boost-invariant longitudinal expansion 

In this section, we study purely longitudinal expansion (in the z-direction) with boost- 
invariant initial conditions |18| . In this case, physics is constant along the space-time hyper- 



bolae t = \Jt 2 — z 2 = const.. Our aim is to establish how the parton equilibration process 
is affected when the full phase-space distribution ([[]) is used instead of the factorized distri- 
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bution (0). We will furthermore compare the time evolution of the parton densities and the 
entropy in the SSPC and pQCD scenarios, and investigate the dependence of equilibration 
on the initial proper time r and the strong coupling constant a s . 



4.1 Factorized vs. full phase-space distribution 

For boost-invariant longitudinal expansion, energy-momentum conservation (|TT|) reads 



e + 



e + p 







(31) 



with e = de/dr, r being the proper time. The rate equations fll~9[ ) - (^1|) assume the form 



n g + — 

T 

■ . n Q 

n a H 



n, 



R 3 n q 1-^ -2R 2 nJl 



n q + — 



R2 n„ 



R2 n„ 



n 9, 
n q rig h 2 g 
n q n q n 2 gi 
n q rig h 2 g 
n„ fin ril 



n q rig n 2 g 
n q hgn 2 g/ 



(32) 
(33) 
(34) 

"q '"q n g) 

M]) are four ordinary differential equations in the variable r, containing 
They are numerically solved with a standard Runge- 
Kutta integration routine. At this point note that, while the evolution of the fluid energy 
and momentum densities are completely decoupled from the evolution of the parton densities 
on account of the equation of state @, energy used up in parton production will be reflected 
by the temperature falling faster than in complete (local) thermodynamical equilibrium. In 
the latter case, T ~ r -1 / 3 



Equations (|3T|) - 
four unknowns, T, X g , X q , and Xg. 



lc l 



For the factorized distribution, the evolution equations fl31|) - fl34|) can be further simpli 



fied 



Xg + b(Xg + Xg) +A T + 4 



Xg + b(X q + X q 



T 3t 



Xg 3T 

A, T 



1 

+ - 

T 



R 3 (1-A 9 )-2R 2 



1 - 



X q Xg 



+ 



Ar, 



T 
3T 



+ - 

T 
1 

+ - 

T 



Ol / X„ X q 



a 1 



X„ A„ 



#2 ^ ^ - ^ 



A, 



(35) 
(36) 
(37) 
(38) 



with 6 = b 2 /a 2 = 21JV//64. 

We are now in a position to check how the equilibration process differs when using 
the factorized phase-space distribution function (@) as compared with the full distribution 
function (0). Fig. ||] shows the proper time evolution of quark and gluon number densities, 
rig and n g , normalized to the corresponding values for A» = 1, h q and n g , for SSPC initial 
conditions at RHIC and LHC (cf. Table []]). The results for the factorized distribution 
function tend to slightly overestimate the degree of equilibration, but the deviation is of the 
order of a few percent only. This was to be expected from Fig. ||. 
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Figure 3: Parton densities for the SSPC RHIC (left) and LHC scenarios (right), normalized 
to the corresponding values for A, = 1 (denoted as hi in the text). The solid line shows 
results for the full phase-space distribution function, while the dashed line those for the 
factorized distribution function / fac . 



13 



1 . o 

. 9 

. 8 

. 7 

O . 6 

. 5 

O . 4 

. 3 

O . 2 

. 1 

. O 

1 . 
O . 9 
. 8 
O . 7 
. 6 
O . 5 
O . 4 
O . 3 
O . 2 
O . 1 
. O 



RHIC : n/n (X = 7 ) 



SSPC 
pQCD 



2.0 3.0 
T [fm/c] 



RHIC 



_/n_CA_= 7 ) 



SSPC 
pQCD 



2.0 3.0 
T [fm/c] 



1 . 

O . 9 

O . 8 

O . 7 

O . 6 

O . 5 

O . 4 

O . 3 

O . 2 

O . 1 

O . O 



I 


I 


i i 


- 


RHIC 


■ n /n ( A = 7 ) 
<r q v q ' 


- 




- 

sspn 






pQCD 




i 





2.0 3.0 
T [fm/c] 




SSPC 
pQCD 



2.0 3.0 
T [fm/c] 



Figure 4: RHIC SSPC and pQCD scenarios (see Table [T]): proper time evolution of the gluon 
(upper left panel), quark (upper right panel), and antiquark densities (lower left panel), and 
the entropy (lower right panel). The densities are normalized to their corresponding values 
for Aj = 1, the entropy (~ st) to its initial value. 



4.2 Comparison of SSPC and pQCD scenarios 

We now compare the proper time evolution of the parton densities and the entropy in the 
SSPC and pQCD scenarios. We exclusively use the full distribution function (TJ) in the 
following. The entropy in a given rapidity unit is dS/drj = sttiR 2 , where ttR 2 is the 
transverse area of the expanding system. Since this area is constant in a purely longitudinal 
expansion, we may consider the product of entropy density and proper time, s r, as a measure 
for the entropy per rapidity unit. 

In Fig. |] we show the parton equilibration process for RHIC initial conditions (see Table 
[l]), for the SSPC model (solid lines) and pQCD (dashed lines). For the SSPC case, the 
hadronization energy density = 1.45 GeVfm -3 is reached at a proper time = 4.15 fm/c, 
at which we stop the time evolution. The relative gluon density, n g /h g , reaches about 0.7, and 
the relative quark density, n q /h q , about 0.5. The entropy increases 13 % before hadroniza- 
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Figure 5: As in Fig. |J for LHC initial conditions. 



tion. For pQCD initial conditions, the QGP phase has little time to evolve, hadronizing 
only about 0.7 fm/c after thermalization. The partons have no time to equilibrate, and the 
entropy increases only ~ 5 %. 

While the SSPC scenario was considered to be net-baryon free, the pQCD scenario shows 
a slight difference in the quark and antiquark initial densities, corresponding to an initial 
baryon number density of about 0.12 fm -3 . We checked that our Runge-Kutta solver respects 
baryon number conservation, which in the purely longitudinally expanding geometry implies 
that the product n# r = {n q — ra 5 )r/3 is constant throughout the expansion. As a consequence 
of nonzero net-baryon number, the equilibrium ratio of n q /n q will be larger than 1, while 
that of riq/hq will be smaller than 1, because n, is computed with Aj = 1 instead of the 
correct equilibrium value A° q = 1/A^ q > 1. For a small initial baryon number density of 
about 0.12 fm~ 3 , however, the deviation of A^ q from 1 is negligible. 

The LHC case is shown in Fig. [|. The parton species are seen to approach equilibrium 
after about 10 fm/c. We note that this time is well before a transverse rarefaction front can 
penetrate to the center of the QGP (see discussion in section |5|). The reason why partons 
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equilibrate is that the system starts at much higher temperature and higher initial values for 
the fugacities than at RHIC (see Table [1]). Particularly noteworthy is that, in the pQCD case, 
the gluons are almost completely equilibrated already at r , while the quarks and antiquarks 
are not. Consequently, quark-antiquark production processes drive the gluons temporarily 
out of equilibrium. Entropy production is of the order of 13% (SSPC) to 17% (pQCD). 



4.3 Sensitivity to To and ot s 

We now investigate the sensitivity of the equilibration process to variations of the initial 
proper time, r , and the strong coupling constant, a s . As discussed above, the values for 
the initial time To in Table [I] are most certainly lower bounds. We therefore consider values 
t > 0.25 fm/c for the SSPC scenario and r > 0.1 fm/c for the pQCD approach. In 
varying r , we keep the produced transverse energy, dE T /drj, and the parton numbers, 
dNi/dr], constant. By eq. (0), the energy density and parton number densities must then 
decrease like ~ 1/tq. In essence, this means that we allow the system to evolve without 
doing longitudinal work. This is certainly an idealization: only when the pressure vanishes, 
the system does not perform any longitudinal work, but in our case pressure builds up as 
the system approaches kinetic equilibrium. 

To facilitate the presentation, we only give the final values for the parton densities and 
the entropy at the end of the expansion of the QGP phase. The evolution of the QGP is 
terminated either at the hadronization time, r^, when the longitudinal expansion has cooled 
the system down to an energy density e/j = 1.45 GeV/fm 3 , or at the rarefaction time, r rarefac , 
when a transverse rarefaction wave reaches the center of the system (see discussion in section 

§■ 

The time evolution equation for the energy density, (|31~1), with the equation of state (H), 
has the solution e/eo = (t /t) 4 / 3 . According to eq. (|30l), the product eoT = (d-Er / 'drj) / '(vri? 2 ) 
is constant for constant transverse energy per unit rapidity and constant transverse area. 
Consequently, the hadronization time Th grows with the initial proper time To according to 
Th = Tg 1 / 4 [(dEx /dn) / (nR 2 €h)] 3 ^ ■ The time spent in the QGP phase is therefore 

_ 1/4 ( dE T /d V \ m 

Ar = r h -r = r (-^-| - r . (39) 

This time increases for small To, has a maximum at Tq = (dEx/ drj) / (A 4 ^T:R 2 €h)i an d then 
decreases again. For the SSPC model at RHIC, Tq ~ 1.67 fm/c, while for pQCD at RHIC, 
Tq ~ 0.27 fm/c. At LHC energies, the values for Tq are quite similar, for the SSPC model, 
t * ~ 11.54 fm/c, and for pQCD, t * ~ 11.48 fm/c. 

The hadronization time grows proportional to the initial transverse energy. For LHC 
energies, the transverse energy is so large that a transverse rarefaction wave (see section 
[5]), travelling with sound velocity c s = l/\/3 into matter at rest, reaches the center of the 
system before the longitudinal expansion has cooled matter down to e^. At z = 0, this 
transverse rarefaction wave reaches the center at time T rarefac = t + R/c s . For a Pb nucleus, 
R = 1.12 A 1 / 3 ~ 6.6 fm and T rarefac ~ t + 11.5 fm/c. While we use r h to terminate the 
time evolution at RHIC energy, for the LHC case, the time evolution is terminated at T rarefac 
instead of Th- 
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Figure 6: Hadronization values of the relative densities and the entropy for a s = 
0.2, 0.3, . . . , 0.6, for the SSPC model with RHIC initial conditions. Upper left: gluons, 
upper right: quarks, lower left: antiquaries, lower right: sr/s r . The dashed line is the time 
the system spends in the QGP phase. 
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Fig. |H shows the case for the SSPC model with RHIC initial conditions. The hadroniza- 
tion values of the densities are closer to equilibrium for larger values of a s . This is obvious, 
since then the right-hand sides of the rate equations ( j32|) - (0) are larger, cf. eq. (|22|) , driv- 
ing the system faster towards equilibrium. Consequently, also entropy production increases 
with a s . 

On the other hand, the hadronization values of the densities are further away from equi- 
librium for increasing values of the initial time To- The reason is that, according to eq. fl30|), 
the values for the initial energy density and the parton densities decrease for increasing r 
(the parton number rapidity density dNi/drj and transverse energy rapidity density dE T /drj 
are kept constant). Surprisingly, this does not have an effect on the initial temperature 
which is to all intents and purposes independent of tq. To understand this, consider the 
factorized expressions (|T0|) and (|17D; reducing energy and parton densities by the same fac- 
tor is achieved by reducing the fugacities only, keeping the temperature constant. However, 
the reduction of the fugacities leads to smaller initial values for the parton densities. This 
puts the initial densities further away from their equilibrium values. This difference is, with 
increasing To, increasingly harder to overcome during the lifetime At, eq. (|39|), of the QGP 
phase. The amount of entropy produced in the RHIC scenarios is proportional to the lifetime 
of the QGP phase, i.e., the time during which chemical reactions drive the system towards 
thermodynamical equilibrium and thus increase the entropy. Note that equilibration is never 
complete at RHIC, unless one uses rather large values of a s . 

The pQCD RHIC scenario is depicted in Fig. |7]. The parton densities behave as in Fig. 
|5|. The main quantitative difference is that, due to the small initial values of the energy and 
parton densities, even for large values of a s the system never reaches chemical equilibrium. 
The behavior of the entropy follows again that of the lifetime of the QGP. 

The LHC case is shown in Figs. [8] and |9|. Despite different initial conditions, the relative 
densities at r rare f ac are remarkably similar in both the SSPC model and pQCD. Again, this 
ratio increases with increasing a s and decreasing r . The entropy increases with a s and 
tq. However, in this case the time the system spends in the QGP phase is constant, At = 
Trarcfac — To = R/c s — 11.5 fm/c for Pb. The increase in entropy can only be explained by 
the fact that the time integral over the right-hand side of eq. (|29|) is larger, if the system is 
further away from equilibrium (which is the case for larger values of tq). As already seen in 
Fig. |5|, entropy production is stronger in the pQCD case. Equilibration is nearly complete 
for large values of a s , independent of the value of tq. The only scenario where the QGP will 
not reach full chemical equilibration is when a s is small, and the initial time is large. 



5 Longitudinal and transverse expansion 

In this section, we consider (cylindrically symmetric) transverse expansion in addition to 
boost-invariant longitudinal expansion. In this case, the time evolution is no longer given 
by ordinary differential equations, and one has to resort to numerical algorithms which solve 
the equations of fluid dynamics, eqs. (|TTD, and the rate equations ([Tj]) - ([Zip. The numerical 
scheme to solve the fluid-dynamical equations used in this analysis is the relativistic Harten- 
Lax-Van Leer-Einfeldt (RHLLE) algorithm [19|, with geometrical corrections performed 



using Sod's method pD|. Since the rate equations are conservation equations with source 
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terms, they can also be solved with the RHLLE algorithm, with the right-hand sides being 
treated with Sod's method. Because of boost-invariance, it is sufficient to solve the equations 
in the central plane, z — 0, where r = t. 

Our results are exclusively for the SSPC model and for the initial values given in Table |IJ 
we did not vary tq = 0.25 fm/c or a s = 0.3. We first show results for a box profile ~ Q(R—r), 
and then for a wounded- nucleon profile ~ 3/2^/l — r 2 /R 2 . To facilitate comparison with 0, 
for the box profile we use R = 7 fm both at RHIC and LHC energies. For the wounded- 
nucleon profile we use R = 1.12 A 1 / 3 to compute the actual values for the radii of Au nuclei 
(at RHIC) and Pb nuclei (at LHC). 



5.1 Box profile 

In Fig. |H] we show the evolution of the system in the t — r plane for RHIC initial conditions. 
In order to highlight the consequences of transverse expansion and facilitate comparison 
with |7], [H, we plot constant energy contours e(r,t) = eo/iV 4 / 3 . If there were no transverse 
expansion of the system, the energy contours for N > 2 would be parallel to the N — 1 
contour. We see that a rarefaction wave travels into the system from the surface, but does not 
reach the center before the interior cools below the hadronization energy density at = 
4.15 fm/c. This agrees with the Runge-Kutta analysis of the purely longitudinal expansion 
presented above. It disagrees with Fig. 1 of |]7j , where the hadronization time is of the order 
of 5 fm/c. Upon hadronization, the relative gluon density, n g /h g , has reached roughly 0.7, 
while the relative quark (and antiquark) density is n q /h q ~ 0.5. Modest transverse velocities 
develop near the surface of the system. We find them to be smaller than those in Fig. 1 
of [[F]. This difference is due to the different numerical scheme used in (the SHASTA 



algorithm of [£Tj; we use the RHLLE algorithm |19|). In contrast to J7J, Fig. 1, we do not 
see that these velocities drive the fluid locally away from chemical equilibrium; in our case, 
the respective fluid elements continue to equilibrate. 

Let us give a simple argument as to why this must be the case. Consider a fluid element 
on a given energy density contour in the t — r plane (upper left panel of Fig. |K]). Then, 
as the system expands, d • u > 0, the energy density drops as a function of proper time, 
u-de = — (e + p)d-u < 0. Since any other contour lying outside the original contour (relative 
to the origin) corresponds to a smaller energy density, it also corresponds to later proper 
times. Therefore, in going from the first contour to the second, the proper time in the rest 
frame of a fluid element increases. Since the actual parton number densities are smaller 
than their equilibrium values, and consequently the right-hand sides of the rate equations 
positive, Ri > 0, we conclude from eq. (|24"D that rii/nf 1 must increase in going from the 
first contour to the second. In order to compare with rii/hi, note that hi is computed at 
the same temperature T as n i; with the fugacity \ set to 1. On the other hand, while Af 1 
is also equal to 1 for the SSPC case considered here (on account of vanishing net-baryon 
density), the temperature T eq in r^ q is in general different from T. Let us assume that 
the initial temperature for an evolution in complete local equilibrium is the same as that 
for the evolution including chemical equilibration. This has the consequence that initially 
hi = nf 1 , and thus rii/hi = rii/nf 1 at r . As noted earlier in section 4.1, in complete local 
equilibrium the temperature falls less rapidly than in chemical non-equilibrium, T eq > T at 
all times r > r . Therefore, rii/hi = (n, / nf 1 ) (T eq / T) 3 > nj/n^ q , with the equality holding at 
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Figure 10: SSPC RHIC scenario: evolution in the t—r plane. Upper left panel: energy density 
contours. The dashed line represents the advance of the rarefaction front. Upper right panel: 
relative density of gluons. Lower left panel: relative density of quarks (antiquarks identical). 
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r . Thus, not only does nj/fij not decrease, it increases even faster than rii/n^ during the 
evolution. Consequently, in the rii/hi — r diagram (upper right and lower left panels of Fig. 
|TO"D , the second contour must lie outside the first contour (with respect to the origin). This 
is exactly what we find numerically. 

Fig. |TT| shows the SSPC LHC scenario. In this case, a rarefaction wave has enough time 
to penetrate into the center of the system, albeit only shortly before the energy density 
falls to th = 1-45 GeV/fm 3 . Until the rarefaction wave reaches the center at proper time 
r = t + R/c s , the evolution of the system at r = is given by a purely longitudinal 
expansion, as discussed in the previous section. Once again, the results obtained with the 
RHLLE algorithm agree well with the analysis using the Runge-Kutta method. As seen 
above, equilibration of the parton species is nearly complete. 

The hadronization time agrees well with that found in [0, Fig. 3. However, in contrast 
to that analysis we again find that the transverse velocities developing in the rarefaction 



23 



14.0 
13.0 
12.0 
11.0 
10.0 

9 . 

8 . 

7 . 

6 . 

5 . 

4 . 

3 . 

2 . 

1 . 



I I 

t [fm/c] 

- 


i i i i i 


- 




e(r,t) = 425 /H ' GeV/fm 










\\\\ 








^^^hadron i zat i on (N-71) 




\\ 
























\V\vv\ 




N=21 


\W\\\ 












\\Vs\A\ 




N= 7 7 


x \\\VV 






















N=7 






i i 


i i i i r 





1 . o 

O . 9 

O . 8 

O . 7 

O . 6 

O . 5 

O . 4 

O . 3 

O . 2 

O . 1 



1 1 1 1 


N~71 








N=21 






N • 1 






N=7 / 




n g/ n g( x g =1 ) along e(r,t) 


= 425/N 4/3 GeV/fm 3 


i i i i 







1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 0.0 1.0 2.0 3.0 4 5.0 6.0 7.0 8.0 
r [fm] r [fm] 



1 . O 

O . 9 

. 8 

. 7 

O . 6 

O . 5 

O . 4 

O . 3 

. 2 

. 1 



~i r~ 



N=37 



N=7 7 




n </ n q ( x q =1 ) along 
e(r.t) = 425/N 4/3 CeV/f, 



_i i_ 



1 . o 

O . 9 

O . 8 

O . 7 

O . 6 

O . 5 

O . 4 

O . 3 

O . 2 

O . 1 



I I I I I 






v r along e(r , t) =425/N 4/3 


GeV/fm 3 «g 




- v r /c 






hadronization (N-71 ) 


















is i i i 


/N= 7 7 





1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 
r [fm] r [fm] 



Figure 11: As in Fig. [10], for the SSPC LHC scenario. 
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front do not impede chemical equilibration. In particular, we do not find the tendency to 
drive the system out of equilibrium at larger transverse distances r, and to overshoot the 
equilibrium values close to the origin, as in Fig. 3 of 0. As explained above, in the Ui/fii — r 
diagram, contours corresponding to smaller energy density (larger proper times) must lie 
outside those corresponding to larger energy density (smaller proper times) relative to the 
origin. 



5.2 Wounded-nucleon profile 

Throughout the analysis presented above a box-profile function has been assumed. While this 
is instructive and facilitates direct comparison between the results from a purely longitudinal 
expansion and those including transverse expansion (as long as the transverse rarefaction 
wave has not yet reached the center of the system), a box profile is clearly an idealization. 

In this section we repeat the analysis of the SSPC RHIC and LHC scenarios with a more 
realistic nuclear profile function, the so-called wounded-nucleon profile 




7Mr) = -Wl--, (40) 



where r is the transverse coordinate. We generate an r-dependence of our initial energy and 
parton density profiles by multiplying the initial values of Table with eq. ((5]). Clearly, 
this implies higher initial densities at the center of the plasma. However, the fluid-dynamical 
expansion begins with a non-zero density gradient already at tq, which, without the factor 
3/2 in (fM]), would lead to faster cooling than a constant nuclear profile function. 



In Fig. [12| we present the analogue of Fig. [10], now computed with an initial wounded- 
nucleon density profile. We observe that the lifetime of the plasma phase is actually larger 
than in Fig. [II], r/, = 4.70 fm/c, instead of = 4.15 fm/c. Due to the initial density 
gradients, transverse velocities develop now over the entire volume of the system. The 
cooling effects of the flow become apparent when considering purely longitudinal expansion 
with an initial energy density of 3/2 x 61.4 GeV/fm 3 . In this case, the QGP lives longer, 
Th = 5.63 fm/c. 

The increase of the initial energy and parton densities by a factor of 3/2 in the center 
of the system has the effect that the partons at small and intermediate r equilibrate a little 
further in the wounded-nucleon scenario than in the box-profile scenario. The relative gluon 
density reaches values around 0.8 (0.7 for the box profile) and the relative quark density 
values around 0.6 (compared with 0.5 in the box profile case). In general, however, the 
box-profile and the wounded-nucleon-profile relative parton density calculations agree to a 
remarkable extent. 



In Fig. [13] we show the LHC case, where the effects of the wounded-nucleon profile 
become more apparent. The plasma cools markedly more rapidly than in the box-profile 
case, hadronizing at a proper time Th = 9.37 fm/c. This is more than 3 fm/c earlier than 
for an initial box profile, in which case Th is essentially moments after r rare f ac , the arrival of 
the rarefaction wave at r = 0. 

While the plasma's lifetime is significantly reduced, parton chemistry, driven by high 
temperatures and high gluon densities, is sufficiently rapid to equilibrate the plasma prior 
to hadronization. Qualitatively, the transverse velocities along contours of constant energy 
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Figure 12: As in Fig. ELO, with an initial wounded-nucleon profile. 
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Figure 13: As in Fig. O, with an initial wounded- nucleon profile. The dotted line in the 
upper left panel depicts the hadronization hypersurface for an initially box-profiled plasma 
(of R = 7 fm), shown in Fig. [TTJ. 
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density in Figs. 12 and 13 agree with the analysis given in ||. No parton evolution plots are 
shown in ||. 



6 Conclusions 

We have investigated chemical equilibration of quarks and gluons at RHIC and LHC energies. 
Assuming that a kinetically equilibrated quark-gluon plasma is created at these energies, we 
have used ideal fluid dynamics coupled to rate equations for the parton densities || 0, ^ to 
study chemical equilibration in the further evolution of the system. Our work is motivated by 
the upcoming experiments at RHIC in the fall of this year, and by recent pQCD estimates 
on parton production in the initial stage of ultrarelativistic nuclear collisions Q. These 
estimates serve as initial conditions to the fluid-dynamical as well as the rate equations. 
Besides pQCD initial conditions, we also considered those given by the so-called SSPC 
model H, in order to compare our results to previous studies. At RHIC, the SSPC model 
constitutes an upper limit for the initial energy and parton densities, while pQCD sets a 
lower limit, as the effect of soft background fields is neglected in this approach. It was 
recently argued |22| that the transverse momentum cutoff po should be chosen to match the 
saturation scale, p sat ~ 1 GeV for RHIC, rather than ~ 2 GeV chosen here. Then, the 
initial conditions obtained from pQCD are rather similar to those in the SSPC approach. 
The produced transverse energy density can also be computed within classical Yang-Mills 
theory p3"l ; the values obtained in this way support the SSPC estimates. 



Our results can be summarized as follows. Chemical equilibration of partons is never 
complete at RHIC energies. The highest degree of equilibration can be reached with SSPC 
initial conditions and an (unphysically) large value of the strong coupling constant, a s ~ 0.6. 
With initial conditions from pQCD, the system never comes close to equilibration. 

The situation changes completely at LHC energies. Here, the initial energy and parton 
densities are already so large that chemical equilibration is complete in almost all scenarios 
considered. Only if a s ~ 0.2 or smaller and the initial proper time r is large, equilibration 
was seen to be incomplete. 

Within the macroscopic transport approach employed here, entropy production due to 
chemical non-equilibrium processes was found to be large, in some cases ~ 30%. This 
has implications for estimating the initial conditions in ultrarelativistic nuclear collisions. 
Commonly, one assumes that the evolution of the system is entropy conserving. Then, since 
s ~ const, x n, the initial entropy ~ SqTq is estimated from the total multiplicity of hadrons 
in the final state by the formula diV^/d^ ~ const, x sqt ttR 2 . However, if entropy grows 
due to chemical reactions, this formula would severely overestimate the entropy in the initial 
state. 

We furthermore established that the use of a factorized phase-space distribution as in 
U |7], [5] is permissible as the error is of the order of a few percent only. We also gave an 
analytical proof, as well as numerical evidence, that transverse flow does not drive the system 
away from chemical equilibrium, in contrast to the results of [/J. 

Future studies should focus on the further evolution of the system into the mixed and 
hadron phases. To this end, previously developed non-equilibrium models could be res- 
urrected |24|. Further applications of the methods presented here include the study of 
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strangeness and charm equilibration at RHIC and LHC energies [12\. Besides this, the re 



action rates should be improved using the full phase-space distribution function instead of 
the factorized one. Other improvements of the expressions for the reaction rates include the 
running of the strong coupling constant. As shown in Ref. | 25| , equilibration is faster when 
the scale of the running coupling constant is allowed to vary with the temperature T. 
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